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AN  ENTROPY  MAXIMUM  PRINCIPLE  AND  RELAXATION  PHENOMENA 


I.  Introduction 

In  this  paper  we  consider  a  general  approach  to  the  description  of  the 
relaxation  of  an  open  system  based  on  a  maximum  Shannon  entropy  principle.* 
The  definition  of  the  Shannon  entropy  involves  an  integral  over  the  full 
temporal  epoch  of  the  relaxation  and  will  be  designated  the  epoch  entropy. 
The  definition  of  the  epoch  entropy  is  based  on  a  dimensionless  quantity 
termed  the  incremental  relaxation  density  (IRD).  The  dimensionless  character 
of  the  IRD  is  maintained  by  the  introduction  of  a  unit  of  time,  and  the  inte¬ 
gral  of  the  IRD  over  the  full  epoch  is  unity  so  it  is  normalized.  The  IRD 
itself  is  defined  in  terms  of  the  approach  to  equilibrium  of  the  singlet 
probability  for  a  particular  state  in  a  discrete  state  space  used  to  describe 
the  system  of  interest.  The  singlet  probability  for  this  state  is  assumed  to 
approach  its  equilibrium  value  monotonically ,  and  at  a  rate  that  is  equal  to 
the  slowest  among  all  states.  When  such  a  state  exists,  it  clearly  controls 
the  ultimate  approach  of  the  system  to  equilibrium,  and  the  relaxation  is  said 
to  be  simple.  We  deal  only  with  simple  relaxation  in  this  paper. 

The  IRD  is  defined  in  such  a  way  that  it  is  positive  with  negative  slope 

throughout  the  epoch  of  relaxation.  It  therefore  has  the  properties  of  a 

2 

Lyapounov  function,  and  provides  a  measure  of  the  status  of  the  relaxation 
process.  Furthermore  the  IRD  is  normalized  so  that  it  also  has  the  mathemati¬ 
cal  properties  of  a  probability  density.  It  is  therefore  an  appropriate 
quantity  to  be  used  in  the  definition  of  a  Shannon  (epoch)  entropy  which 
provides  a  measure  of  information  about  the  relaxation  process  throughout  the 
epoch  for  the  system  of  interest. 

We  maximize  the  epoch  entropy  subject  to  minimal  constraints  for  a  relax¬ 
ation  process,  namely,  that  the  IRD  is  normalized,  and  that  there  is  a  mean 
Manuacript  approved  April  30,  1984. 
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time  of  relaxation  that  corresponds  to  the  first  moment  of  the  IRD.  The 
resulting  form  for  the  IRD  that  assures  maximum  epoch  entropy  is  a  linear 
exponential.  Since  it  is  derived  directly  from  a  maximum  entropy  principle 
subject  only  to  the  minimum  constraints  necessary  to  specify  a  relaxation 
process,  we  consider  the  linear  exponential  to  be  the  fundamental  form  for  the 
IRD.  With  this  IRD  in  hand,  we  can  find  a  general  expression  for  maximum 
epoch  entropy.  However,  we  note  that  the  time  scale  on  which  the  IRD  takes  a 
linear  exponential  form  is  not  in  general  the  time  scale  on  which  measurements 
are  made.  We  argue  that  the  value  of  maximum  epoch  entropy  should  not  depend 
on  the  time  scale  that  is  used  to  describe  it.  We  are  led  thereby  to  enunci¬ 
ate  a  principle  of  invariance  of  maximum  epoch  entropy  in  which  the  maximum 
epoch  entropy  evaluated  in  terms  of  the  transformed  IRD  is  equal  to  the  maxi¬ 
mum  epoch  entropy  evaluated  in  terms  of  the  fundamental  IRD.  This  principle 
leads  to  an  invariance  relation  between  units  of  time  on  the  fundamental  time 
scale  and  a  transformed  time  scale. 

We  then  use  some  general  arguments  to  rationalize  the  use  of  a  monomial 

(with  positive  fractional  power)  time  scale  transformation  in  the  description 

of  relaxation  phenomena  in  an  open  system.  The  fundamental  form  corresponds 

to  the  case  when  the  monomial  is  at  the  linear  limit.  However  the  general 

3 

monomial  transformation  leads  to  an  empirically  well-established  description 

of  relaxation  phenomena  that  also  coincides  with  a  description  based  on  an 

4  5 

existing  microscopic  model  of  relaxation.  ’  Furthermore  the  invariance 
relation  between  units  of  time  that  is  derived  from  the  principle  of  invari¬ 
ance  of  maximum  epoch  entropy  leads  to  renormalization  relations  for  scaling 
parameters  and  activation  energies  that  enter  into  the  physical  manifestations 
of  these  units.  Similar  renormalization  relations  also  arise  in  the  aforemen- 

and  have  been  repeatedly  verified  empirically.^ 


tioned  microscopic  model, 


II.  The  Incremental  Relaxation  Density  (IRD)  and 

Maximum  Temporal  Epoch  Entropy 

We  considered  a  perturbed  open  system  that  is  relaxing  to  equilibrium 
with  order  parameter  (time)  0=0-0q,  where  0Q  denotes  the  time  when  the  relaxa¬ 
tion  begins.  The  system  is  taken  to  be  characterized  by  a  singlet  probability 

P  (0)  over  a  discrete  state  space  with  states  labelled  by  s.  P  (»)  are  the 
s  s 

equilibrium  values  of  Pg(0)  so  a  state  retains  its  identity  throughout  the 
relaxation  process.  Since  Pg (0)  are  state  probabilities,  we  have  for  all  0 

Fs(e)  =  1 

Among  the  states  s,  there  is  a  subset  consisting  of  at  least  two  states 
for  which  the  relaxation  is  slowest.  These  will  be  called  the  slowest  relax¬ 
ing  states.  The  evolution  of  the  probabilities  associated  with  such  states 
therefore  control  the  ultimate  approach  to  equilibrium.  We  say  that  a  relaxa¬ 
tion  is  simple  if  for  at  least  one  of  the  slowest  relaxing  states,  the  state 
probability  approaches  its  equilibrium  value  monotonically .  Let  us  designate 
that  state  (or  a  chosen  one  of  a  number  of  such  states)  to  be  the  rth  state 
with  probability  Pf(0)  and  equilibrium  probability  Pr(®).  We  now  make  a  weak 
assumption  (to  be  strengthened  below)  that  because  of  the  nature  of  a  relaxa¬ 
tion  process,  Pr(0)  approaches  Pr(°°)  sufficiently  rapidly  that  the  integral  of 
Pr(Q)”Pr(°°)  over  the  full  epoch  of  the  relaxation  process  is  finite. 

Nr  =  J^(Pr(e)-Pr(®)]d0/TQ  (2) 

Here  has  the  unchanging  sign  of  the  difference  Pr(0)-Pr(00) ,  and  Tg  is  a 
(for  the  moment  arbitrary)  constant  unit  of  time  on  the  0-scale  so  that  the 
integral  is  dimensionless. 

We  now  observe  that  for  simple  relaxation,  we  may  define  a  dimensionless 
incremental  relaxation  density  (IRD)  to  be 

f (e)  =  N‘1[Pr(e)-Pr(“)]  (3) 
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The  IRD  is  a  positive  semi-definite  function  of  6  with  negative  slope  that  has 
its  maximum  at  0=0  and  which  vanishes  as  0-*».  Therefore  by  its  positivity  and 
negative  slope  for  0^0,  the  IRD  has  the  properties  of  a  Lyapounov  function. 
In  addition  the  integral  of  the  IRD  over  the  epoch  is  just  unity.  Thus  the 
IRD  also  has  the  properties  of  a  dimensionless  probability  density.  Its 
accumulation  integral 

F(0)  =  ;®f(0,)d6,/x0,  F(o0  =  1  (4) 

is  positive  with  positive  slope  for  all  0^0.  Now  the  properties  of  F(0)  could 
also  arise  from  an  epoch  integral  over  a  density  function  which  is  positive 
but  with  a  slope  whose  sign  changes  during  the  epoch.  It  is  interesting  to 
speculate  that  the  conditions  on  a  Lyapounov  function  may  be  similarly 
weakened  so  that  the  requirement  of  negative  slope  is  replaced  by  a  condition 
that  the  accumulation  integral  of  a  generalized  Lyapounov  function  over  the 
epoch  is  finite.  Of  course,  then  the  relaxation  would  not  be  simple  so  we 
leave  this  point  for  further  consideration  elsewhere. 

We  have  already  observed  that  the  slowest  relaxing  states  and  hence  also 
the  IRD  describe  the  status  of  the  ultimate  relaxation  of  the  system  to  equi¬ 
librium.  Further  as  also  previously  noted,  the  IRD  has  the  mathematical 
properties  of  a  dimensionless  probability  density.  Thus  we  are  led  in  analogy 
with  procedures  used  by  Shannon*  to  define  a  Shannon  temporal  epoch  entropy  S0 
as  a  measure  of  information  about  the  relaxation  process  over  its  full  epoch. 

S0  =  -j^f(0)£n[f(0)ld0/T0  (5) 

It  should  be  noted  that  S0  is  an  implicit  function  of  0q. 

We  now  require  that  there  is  a  minimum  of  information  about  the  relaxa¬ 
tion  process  subject  only  to  those  constraints  appropriate  for  specifying  the 
process  to  be  a  relaxation  process.  Equivalently,  this  means  the  epoch 
entropy  is  a  maximum  subject  to  appropriate  constraints.  We  already  have  the 
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condition  that  the  IRD  is  normalized,  by  its  definition  Eq.  (3).  In  addition, 
a  relaxation  process  is  typically  characterized  by  a  relaxation  time.  We 
specify  the  relaxation  time  to  be  represented  by  the  first  moment  of  the  IRD, 
namely 

<0>  =  J^0f(0)d0/t0  (6) 

This  provides  a  somewhat  more  stringent  condition  on  the  rapidity  with  which 
P^(0)  approaches  P  («)  than  required  for  Eq.  2. 

Now  introducing  two  Lagrange  multipliers  A^  and  in  the  usual  way,  we 
obtain  from  the  variation  of  S0  that  the  maximum  epoch  entropy  occurs  for 


f^(0)  =  expf-CAj  +  A20)]  (7) 

Consistency  with  Eqs.  (3)  and  (6)  allows  us  to  determine  A^  and 

exp(-A^)  =  T0/<0>  ,  A2  =  <0>  1  (8) 

Then  the  IRD  consistent  with  the  maximum  epoch  entropy  has  the  form 

fM(6)  =  [tQ/<0>]  exp [ - (0/ <0> ) ]  (9) 

The  maximum  epoch  entropy  can  itself  be  determined  to  be 

S0  =  1  +  £n<0>-£ntQ  (10) 


M 

so  the  specific  value  of  S0  is  determined  modulo  a  choice  for  T0. 

There  are  several  important  observations  that  should  be  made  concerning 
Eqs.  (9)  and  (10).  First  from  Eq.  (9)  we  note  that  the  linear  exponential 
form  is  consistent  with  the  Lyapounov  function  properties  of  the  IRD.  Also 
the  linear  exponential  form  for  the  IRD  represents  a  fundamental  form  in  the 
sense  that  it  arises  from  the  maximum  epoch  entropy  principle  subject  only  to 
the  minimum  constraints  necessary  to  characterize  the  (simple)  relaxation 
process.  However,  it  has  long  been  recognized^  that  an  exponential  form  is 
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not  completely  appropriate  for  the  description  of  the  dynamics  of  decay  at  the 
shortest  times  (or  at  the  highest  frequencies).  Therefore  the  point  6=0  is 
physically  singular.  Here,  we  understand  this  to  mean  that  there  must  be  a 
pre-relaxation  regime,  before  simple  relaxation  begins,  that  may  be  character¬ 
ized  by  a  high  frequency  cutoff  u)^.  Thus  the  present  discussion  is  valid  only 

for  long  times  corresponding  to  frequencies  lower  than  cu  .  Empirical  evi- 

3  8  1 1 

dence  indicates  that  for  many  cases  ujc  may  be  in  the  range  10-10  Hz  for 

relaxation  phenomena  in  condensed  matter. 

M 

It  is  also  useful  to  note  that  although  f  (0)d0/Tg  is  independent  of  the 
M 

time  unit,  f  (9)  alone  takes  on  a  particularly  simple  form  when  tg  is  taken  to 

M 

equal  to  <9>.  Namely  f  (0)  is  then  equal  to  unity  at  6=0  and  homogeneous  in 

0/<0>.  Thus  <0>  is  a  natural  time  unit  for  the  0  time  scale.  A  similar 

M 

simplification  occurs  for  Sg  which  just  becomes  unity  when  Tg=<0>.  However 

M 

for  any  particular  choice  of  Tg,  Sg  is  an  epoch  integral  whose  value  should  be 
independent  of  the  time  scale  on  which  it  is  described.  This  condition  of 
invariance  serves  to  place  a  stringent  condition  on  the  relationship  between 
tg  and  the  unit  of  time  on  the  new  time  scale.  This  is  probably  obvious  for  a 
linear  time  scale  transformation  but  we  are  interested  here  in  a  wider  class 
of  time  scale  transformations.  We  are  motivated  here  by  the  fact  that  there 
is  no  reason  to  believe  that  the  time  scales  on  which  measurements  are  made 
coincide  with  the  time  scale  on  which  the  IRD,  f  (9),  appears  in  its  funda¬ 
mental  linear  exponential  form. 

III.  An  Invariance  Relation  for  Time  Scale  Units 

A  time  or  order  parameter  is  just  a  positive  cumulative  function  that 
increases  monotonically  from  an  origin.  We  have  already  introduced  the  0  time 
scale,  its  corresponding  difference  time  scale  0=0-0^  and  its  constant  unit 


i 


i 


» 


wm 

i 


« 


of  time  Tq.  In  a  similar  way,  we  introduce  another  time  scale  t  with  differ¬ 
ence  time  scale  t=t-EQ,  and  constant  unit  of  time  From  the  context  of 

application,  it  is  clear  that  we  should  align  the  respective  origins  0q  and  £q 
since  both  represent  the  time  when  the  (simple)  relaxation  process  begins 
(after  a  pre-relaxation  time  regime).  Thus  if 

0=6(t)  ,  0(t=O)  =  0  .  (11) 

We  specify  further  that  d0/dt  is  positive  and  finite  everywhere  except  pos¬ 
sibly  at  isolated  points.  For  example,  such  an  isolated  point  might  be 
expected  at  the  end  of  the  pre-relaxation  regime  which  corresponds  to  the 
origin  for  the  time  regime  of  simple  relaxation. 

We  can  now  describe  the  incremental  relaxation  de  ^  (IRD)  and  maximum 

M 

epoch  entropy  in  terms  of  t.  Namely,  we  obtain  a  tra  turmed  IRD  g  (t)  by  a 

M  M 

direct  transformation  from  f  (0)d0/tg  to  g  (t)dt/xt.  W  find  using  Eq.  (9) 

gM(t)  =  exp{-[0(t)/<0>]Htt/<0>]  [d0(t)/dt]  .  (12) 

The  maximum  epoch  entropy  is  now  taken  to  be  form-invariant  so  that  it  is 

M 

written  in  terms  of  g  (t)  and  an  integral  over  the  total  epoch  on  the  t-scale. 
St  =  J^8M(t)£n8M(t)dtAt  =  l+2n(<0>) 

'  -Cexp["  4>  Un(rt  5t  }  (13) 

Within  the  integral  on  the  right  hand  side,  d0/dt  is  expressed  as  a  function 

of  <0>  and  0/<0>.  We  now  use  the  requirement  that  S'  is  invariant  independent 

of  the  time  scale  used  in  its  evaluation  together  with  Eqs .  (10)  and  (13)  to 

obtain  a  general  invariance  condition. 

2n(tQ/xt)  =  J^exp[-(0/<0>) ]£n[d0/dt}d0/<0>  (14) 


•  ( 
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To  understand  the  ramifications  of  this  relationship  we  should  consider  a 
specific  time  scale  transformation.  The  immediate  problem  then  is  to  choose  a 
time  scale  that  is  meaningful  for  measurements  of  relaxation  phenomena. 

IV.  A  Measurement  Time  Scale  for  Relaxation 

We  observe  that  even  after  a  time  scale  transformation,  exp [ -0(t)/<0> ]  is 

positive  for  all  6(t),  t>0,  has  value  unity  at  0(0)  because  of  Eq.  (11),  and 

vanishes  as  0(t)-»0(»).  It  was  previously  argued2  that  such  a  decay  in  an  open 

2 

system  can  be  described  in  terms  of  the  squared  magnitude  | c ( t ) I  of  the 
antocorrelation  amplitude  of  a  decaying  state  |R>.  That  argument  may  be 
sharpened  as  follows.  A  density  matrix  representing  the  system  at  any  time  is 
a  Hermitian  trace-class  operator  with  eigenvalues  that  are  positive,  real,  and 
nonvanishing,  and  which  lie  between  0  and  1.  The  eigenvalues  belong  to  the 
complete  set  of  the  eigenstates  of  the  density  matrix  which  includes  both 
|  non-stationary  decaying  and  stationary  states.  In  equilibrium  the  decaying 

states  no  longer  have  non-vanishing  eigenvalues  so  that  the  residual  eigen¬ 
values  of  the  density  matrix  belong  only  to  its  stationary  eigenstates.  The 
|  decaying  states  are  therefore  not  relevant  for  the  description  of  equilibrium. 

Hence  the  approach  of  the  incremental  relaxation  density  (IRD)  to  zero  matches 
the  vanishing  of  the  autocorrelation  of  a  decaying  state.  In  such  a  case 

exp [ -0 ( t ) / <0> ]  =  |c(t)l2  (15) 

where 

,  c(t)  =  <R  |exp(- ( iDt ) j R>  ,  t  ^  0  (16) 

The  state  i  R>  is  associated  with  the  continuous  spectrum  of  the  development 
operator  D  and  is  orthogonal  to  all  the  stationary  eigen-states  of  the  density 
*  matrix  operator.  It  is  assumed  that  the  development  operator  has  a  non- 


singular  continuous  spectrum. 


If  pc  denotes  the  spectral  projection  of  D 
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(17) 


D  =  Je  dp£  =  /elexelde  , 

then  the  function  <Rl|J  |R>  is  absolutely  continuous.  Its  derivative 

o 

p(fi)  =  ^  <RIM£IR>  =  <R|£><£|R>  (18) 

can  be  interpreted  as  the  spectral  distribution  of  the  state  IR>.  In  other 

E+AE 

words,  the  integral  Jg  p(e)de  is  the  probability  that  the  eigenvalue  of  the 
state  I  R>  lies  in  the  interval  (E,E+AE).  The  function  p(e)  has  the  following 
properties . 

i)  p(£)  ^  0 

ii)  Jp(e)d£  =  1,  since  <R|R>  =  1  (19) 

iii)  p(e)  =0  ,  £  t  P£ 

Thus  p(fi)  is  a  probability  density.  Now  it  follows  that  c(t)=<R| exp(-iDt)l R> , 
t  ^  0,  can  be  expressed  as 

c(t)  =  /exp(-iet)  p(e)d£  (20) 

Thus  c(t)  is  a  characteristic  function  for  the  spectral  probability  density  of 
the  decaying  state. 

We  now  note  that  there  is  some  arbitrariness  in  the  specification  of  I  R> 
and  its  spectrum.  In  general  we  consider  I  R>  to  represent  a  decaying  state 
for  a  complex  many-body  system.  Then  the  description  of  the  decay  should  not 
depend  on  an  exact  specification  of  the  number  of  components  in  the  system. 
In  other  words  we  should  obtain  the  same  form  of  decay  respectively  for  one 
portion,  or  several  portions  of  the  system.  Mathematically  this  means  that 

p(£)  must  be  superposable  probability  density,  as  discussed  by  Rajagopal  et 

8  9 

al.  so  that  p(£)  must  be  a  stable  distribution.  Recently  Weron  et  al.  have 

considered  the  application  of  the  mathematical  theory  of  stable  distributions 

for  systems  with  semibounded  spectra  to  problems  of  interest  in  physics.  Of 


pertinence  to  present  considerations  is  the  result  that  for  a  semibounded 
stable  distribution,  the  appropriate  characteristic  function  for  causal  time, 
t^O,  has  squared  magnitude 

|c(t)|^  =  exp(-at^),  a>0,  0<b<l,  t^O  .  (21) 

A  comparison  with  Eq.  (15)  then  provides  a  time  scale  transformation  appropri¬ 
ate  for  the  description  of  relaxation  phenomena.  Namely  the  fundamental  time 
scale  is  related  to  the  time  scale  of  measurement  by 

9  =  atb,  a>0 ,  0<b<l,  UO  (22) 

where  a  =  a<0>.  We  shall  apply  this  time  scale  transformation  in  a  discussion 
of  the  invariance  relation,  Eq.  (14),  in  the  next  section.  Before  that, 
however,  some  remarks  concerning  Eq.  (21)  are  appropriate. 

We  note  first  the  form  of  |c(t)l  in  Eq.  (21)  is  consistent  with  the 
Paley-Wiener  theorem  in  fourier  transform  theory. ^ The  Paley-Wiener 
theorem  states  that  a  necessary  and  sufficient  condition  for  the  existence  of 
a  semibounded  p(e)  (with  properties  given  in  Eq.  (19))  which  is  the  fourier 
transform  of  a  square  integrable  c(t)  is  that 

Jl  <!-i&)l  «  <  -  (23) 

2 

This  condition  is  guaranteed  if  for  t-*00,  | c ( t ) |  is  greater  than  or  equal  to 
the  right  hand  side  of  Eq.  (21).  The  important  point  is  that  for  the  stable 
distribution  argument  used  here,  Eq.  (21)  is  an  equality  and  valid  throughout 
the  full  epoch. 

It  should  also  be  noted  that  although  the  limit  b=l  can  never  be  reached, 
b  may  approach  arbitrarily  close  to  1 .  It  is  that  latter  limit  that  corre¬ 
sponds  to  the  fundamental  form.  The  deviation  of  b  from  unity  is  a  manifesta¬ 
tion  of  the  complex  nature  of  the  system  as  it  appears  on  the  time  scale  of 


measurement. 


V. 


Parameter  Renormalization 


We  now  consider  the  effect  of  the  time  scale  transformation  Eq.  (22)  as 
applied  in  Eq.  (14).  We  obtain 

M  ^  ^  [<e>/a)1/b}  =  -(1  -  I)Y  (24) 

t 

where  y  =  0.577215  ...  is  the  Euler  constant.  If  we  now  choose  Xg  to  coin¬ 
cide  with  the  natural  unit  on  the  0-scale,  namely  <0>,  we  obtain  the  relation¬ 
ship 

<0>  =  abbexp[  (l-b)y]xb  (25) 

It  is  convenient  to  define  a  new  time  unit  on  the  t-scale,  namely 
Tt  =  b  exp [  (l-b)y/b]tt  (26) 

so  that 

<0>  =  a  x°t  (27) 

Then  using  Eq.  (12),  we  obtain  for  g^(t)dt/xt 

gM(t)dt/xt  =  exp(-(t/xt)b]  b(t/xt)b  1  dt/xt 

=  iM(t)dtdt  (28) 

The  function  g  (t)  is  just  the  incremental  relaxation  density  (IRD)  on  the 
t-scale  with  time  unit  xt> 

It  is  useful  to  relate  this  IRD  to  a  relaxation  rate  W(t)  to  facilitate 
comparison  with  the  usual  description  of  relaxation  phenomena.  Note  that  the 
unrelaxed  portion  at  any  time  t  is  just 

1  “  Jb8^(t)dt/it  =  exp[-(t/x)b]  (29) 

-M 

The  relaxation  rate  in  a  particular  time  interval  is  just  the  ratio  of  g  (t)/xt 
to  the  remaining  unrelaxed  portion,  namely 
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W(t)  =  b(t/Tt)b_1(Tt)"1  ,  0<b<l  (30) 

We  have  previously  discussed**  such  time  dependent  relaxation  rates,  and  used 
a  phenomenlogical  introduction  of  a  monomial  time  scale  transformation  for 
their  application  in  the  description  of  relaxation  phenomena.  On  the  t-scale, 
the  relaxation  rate  decreases  with  time.  Hence  for  measurements  made  on  a 
t-scale,  the  relaxation  will  appear  to  have  a  long  tail. 

The  decay  represented  by  Eq.  (21)  and  the  relaxation  rate  in  Eq.  (30)  are 
the  prototypical  forms  observed  in  many  types  of  relaxation  experiments  in  a 
wide  range  of  condensed  matter  samples.  Even  more  remarkable  are  the  applica¬ 
tions  of  Eq.  (27).  Its  significance  resides  in  the  fact  that  the  t-scale  is 
only  the  time  scale  of  measurement  and  not  the  fundamental  time  scale  for  the 
maximization  of  epoch  entropy  subject  to  the  minimum  constraints  necessary  to 
specify  a  simple  relaxation  process.  Thus  parameter  dependences  determined 
for  It  are  only  apparent.  The  actual  dependencies  are  those  that  enter  into 
<0>.  The  two  dependencies  are  related  by  Eq.  (27). 

For  example,  consider  a  material  made  of  polymer  with  a  relaxation  time 
dependent  on  the  molecular  weight,  M,  of  its  chemical  building  blocks.  On  the 
fundamental  scale,  we  take  the  relaxation  time,  for  specificity,  to  have  a 
simple  monomial  dependence  on  M. 

<0>  =  kf#  (31) 

Then  on  the  time  scale  of  measurement,  the  relaxation  time  tt  has  the  depen¬ 
dence 

Tt  «  MP/b  (32) 

T  therefore  has  a  power  dependence  that  appears  greater  than  it  actually  is 
since  0<b<l.  If  now  we  consider  polymer  samples  with  different  values  of  M, 
relaxation  measurements  will  typically  provide  different  values  of  the  frac- 
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tional  exponent  b.  The  prediction  of  the  invariance  relation  is  that  the 
M-dependence  of  will  vary  in  such  a  way  that  p  is  the  same  for  the  samples 
with  different  M.  On  the  other  hand,  the  measured  (5/b  would  appear  anomalous 
without  an  understanding  of  the  need  for  renormalization.  Similarly  it  may 
also  be  possible  to  change  the  fractional  exponent  b  by  changing  the  environ¬ 
ment  in  which  relaxation  measurements  are  made.  Then  the  prediction  is  that 
the  power  of  M  that  is  measured  for  x  changes  with  b  in  such  a  way  that  p 
again  remains  constant. 

A  more  general  application  occurs  in  the  case  of  the  variation  of  <0>  and 
T  with  temperature.  In  the  usual  way  we  assume  an  Arrhenius  dependence  so 
that 

<0>  =  <6^  exp(EA/kT)  (33) 

Here  is  an  activation  energy  and  <0M>  is  the  relaxation  time  on  the  funda¬ 
mental  scale  at  (nominally)  infinite  temperature.  For  simplicity  sufficient 
for  our  present  expository  purposes,  we  consider  that  the  fractional  exponent 
b  to  be  independent  of  temperature.  On  the  time  scale  of  measurement,  we  find 

Tt  =  (lt)0>exp(E|/kT)  (34) 

Here  E^  is  an  effective  activation  energy.  It  is  clear  that  for  the  invari¬ 
ance  relation  Eq.  (27)  to  hold,  E*  must  be  related  to  the  actual  activation 
energy  EA  by  the  renormalization  relation 

E*  =  EA/b  (35) 

Thus  measured  activation  energies  invariably  seem  to  be  greater  than  they 
actually  are. 

The  remarkable  feature  of  results  like  Eqs.  (21),  (30),  (32)  and  (35)  is 
that  they  and  their  mutual  dependencies  are  consistently  verified  by  experi- 


ment  on  many  kinds  of  condensed  matter  for  a  wide  range  of  relaxation  phenom- 
3 

ena.  This  verification  lends  credibility  to  the  ideas  that  have  been  invoked 
in  the  present  development.  These  include:  the  concept  of  incremental  relax¬ 
ation  density  for  simple  relaxation;  the  principle  of  maximum  epoch  entropy; 
the  existence  of  a  fundamental  time  scale  for  simple  relaxation;  the  need  for 
a  time  scale  transformation  to  obtain  the  time  scale  of  measurement;  the 
invariance  of  the  value  of  the  maximum  epoch  entropy  when  evaluated  on  dif¬ 
ferent  time  scales;  and  the  procedure  presented  in  Sec.  IV  to  justify  the 
choice  of  time  scale  to  correspond  to  the  time  scale  of  measurement  for  relax¬ 
ation  phenomena. 
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